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modified  program  to  treat  steep  fronted  wave  instabilities  in  low  smoke  tactical  motors 
was  verified  by  solving  a  number  of  motor  test  cases  with  disturbance  amplitudes  as  high 
as  80'^  of  the  mean  pressure. 
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The  objective  of  the  research  work  presented  in  this  reixirt  w'as  to  evaluate  existing 
finite  difference  methods  and  select  and  apply  the  one  most  appropriate  for  coni(>utinK  the 
propagation  of  multiple  shock  w'aves  in  a  two-phase  rocket  combustion  chamber. 

In  order  to  be  acceptable  for  the  intended  application,  a  finite  difference  integration 
scheme  must;  preserve  the  high  frequency  content  of  the  wraveforms,  be  relatively 
nondissipative  and  nondispersive  after  many  wave  cycles;  be  capable  of  describing  a  shock 
wave  as  a  sharp  discontinuity;  and  be  capable  of  properly  treating  the  reflection  of  shock 
waves  from  boundaries  and  the  partial  reflection  and  transmission  ot  discontinuities. 

The  finite  difference  schemes  of  the  Split  Coefficient  method,  the  X -scheine,  Kubin 
and  Ifurstein,  MacCormack,  Lax- Wendroff,  Godunov,  Kusanov,  the  Flux-Corrected- I  ransport 
scheme  of  boris  and  book,  Chorin's  implementation  of  Glimm's  Method,  the  Hybrid  scheme 
of  Ffarten  and  Zwas  and  the  Artificial -Compression  method  of  Marten  were  evaluated.  All 
of  these  schemes  are  finite  difference  approximations  to  the  derivatives  arising  in  thr; 
conservation  laws  and  can  treat  an  arbitrary  system  of  conservation  laws. 

As  a  first  step  in  evaluating  the  |x>tential  of  these  schemes,  they  were  utilized  in  the 
solution  of  the  shock  tube  problem  for  the  rxie- dimensional  tiilerian  form  of  the  gas 
dynamic  conservation  equations  for  an  inviscid,  non  heat-conducting  fluid  and  for  the; 
solution  of  the  linear  wave  equation  problem.  The  final  test  problem  was  thr?  solution  of 
the  one-dimensional  nonlinear  hyperbolic  equations  describing  finite  amplitude  wave  and 
shock  propagation  in  a  closed  end  tube.  In  addition,  the  ability  to  S[)ectrally  analyze  the 
computed  results  was  developed.  This  capability  simplifies  the  interpretation  of  the 
coin()lex  waveforms  and  facilitates  comparisons  amrxig  the  various  finite  difference 
integration  schemes. 

The  results  of  ap()lying  several  r)f  the  alxjve  rnentioneo  schemes  tf)  the  last  test  case 
are  presented  and  discussed  in  tliis  report.  It  had  l)een  ctjncludetJ  that  for  tlie  presi-nt 
problem,  a  method  based  up<jn  the  combination  of  the  Lax-Wendrof f.  Hybrid  and  Artificial 
Compression  scheme  was  found  to  bo  superior  to  the  <ither  schemes  tested.  I  his  sc  tienie 
was  incorporated  into  the  nonlinear  instability  program  developcMi  by  I  eviiu'  and  (Tdick  that 
described  the  comlxjstion  and  flow  inside  a  solid  rocki>t  motor.  Finally,  the  ability  (»f  this 
scheme  to  treat  various  initial  disturbanc<;s  in  a  solid  rocket  motor  suen  as  a  first  m(xi<' 
disturbance  and  standing  or  traveling  (tulses,  is  demonstrateii. 
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SLCIION  1 


intkijix^k;  1  i(JN 

lliis  report  presents  the  results  ot  invest l^at ion  desifjneu  to  select  a  satistattor\ 
inethoti  tor  conrputin^t  the  pro()aitat  ion  ot  stecfi-  f  routed  shocklike  tvaveiortos  id  a  tvvo-/;l)asi‘ 
rocket  combustion  chanitx'r. 

Idctical  solid  rocket  motors  are  frequently  subject  to  a  combustion  instability 
problem  at  some  |ioint  in  the  design  cycle.  vdiim  instability  is  encountered,  it  can  take 
one  ot  several  forms,  e.yt.,  linear  or  nonlinear,  lonj^itudinal  or  tangential.  Over  tiie  last 
twenty  years,  considerable  resources  have  been  expended  to  understand,  (ireciict,  control, 
and  eliminate  combustion  instability  in  solid  rocket  motors.  Most  ot  tliis  effort  has  (x'en 
devoted  to  linear  instability  problems,  and  as  a  result,  sucli  proble  ns  can  now  lie  treateii 
in  a  rational,  cost  effective  manner.  In  comparison,  little  work  has  bec^n  acconiplished 
towards  the  understanding  and  resolution  ot  nonlinear  combustion  instability  (irobh'ms. 
Ihus,  when  nonlinear  instabilities  are  encounteriKl,  the  solution  is  tix)  oftt'ii  an  ex/iensivt* 
cut -and -t  ry  [irocess. 

Linear  instafiility  is  clia  racte  rized'  by  small  dm()litude,  simisoidal  oscillat  rrxis  that 
originate  from  the  amplification  of  irrf initt’simal  ranikj/n  disturbances  in  tlie  miotor  cliamlHO. 
Nonlinear  instability  is  usu<illv  characterized  by  larj’e  arriplitude  oscillations  havinp 

steep- f  routed,  shocklike  waveforiris,  and  is  initiated  by  random  linit(“  >imi)litude  I'vents  siirh 
as  the  exjiulsicjn  ot  an  i;;niter  or  insulation  I  raj;ment  through  the  noz/le.  Nonlinear 

instalii  I  it  ies  are  niodeled  usiny.  Ixith  "exact"  and  "approximate"  mathematiral  techniques. 
I  lie  "exact"  methoiis  of  Levine  and  bulick,^  .ind  tiu/se  of  Ktxiker  and  /inn,'^  sei’k  to  solve 
lurnerically  tiie  nonlinear  partial  ditlerential  equ.dions  y,overnm!t  Ixith  the  me.m  and  tinie 

dependent  flow  in  t!u‘  combustion  chamber,  as  well  as  thr-  <<»mbustion  res)X)nsc  ot  the  solid 

'  4 

proix'llant.  I  he  "aptiroximate"  methorts  ot  <  uliiK,  ano  ot  Itiwx'll,  et..d.,  utilize 


I.  I  evirie,  |.  N,  and  l  ulick,  1.  1.  (  .,  "Nonlinear  Xnalysis  ot  Solui  koi  ket  i  omliust  ion 
Instability,"  AIKM  U'chnictil  Kepirt  IIL-JA-ah,  tictolHO 

Z.  Kooker,  I).  I.  aiui  /inn,  b.  I.,  "NumerKal  Solution  ot  \xial  Instabilities  in  SoluJ 
I’ropellant  Kocket  Motors,"  Kith  |  \NN  M  vaimbiistion  Meet  my;,  \ol.  1,  Naval  war  i  olleye, 
Newjxirt,  K.I.,  August  lb7.i,  t  I’lA  I’ub.  Z-t  t. 

t.  t  iilick,  I.  1.  "Nonlinear  lU'havior  ot  Acoustir  Wavr's  in  (omlaistion  i  hampers,’  IDili 

)\NN  \t  (ompiistion  Meeting,  Vol.  I,  Naval  war  <  olli'ge,  Newjxirt,  K.I.,  Aiignsl  l‘)7t. 

4.  Iiowell,  I.  A.,  I'admanaooan,  M.  s.,  arul  /inn,  It.  I.,  "Approximate  Noiilmeai  UiaKsis  ot 

Solid  Kill  ket  Motors  and  I  -  I 'm  me  r  s, "  MKM  -  I  A  -  7  /  -  4<t,  |uK  l't77. 


perturbation  techiii'iufs  anil  tiarnumK  analysis  to  reilui  »■  tne  i;ovi-rninp  ill  1 1 1- riait  la  I 
equations.  I  ,.i  ii  ot  tiiese  netaods  has  certain  advantages,  disailvant ages,  .md  limitations. 
I  he  two  classes  ot  approaches  comp  lenient  each  other,  and  etiorts  to  develop  Imth  Uirther 


are  warranteu. 

1  tie  existing  "exact"  nonlinear  instaliility  programs  ivx're  ilevelopeil  alxiut  seven  vc.irs 
ago  and  are  not  capable  ot  treating  the  multiple  shock  steep- f  ronted  t\|K'  r,t  ins!  aij  tl  it  les 
that  occur  in  reduced  and  ininimum  smoke  tactical  motors  develo|)id  sirue  then.  Uie 
oo/ective  ot  tlie  /iresent  riysearcb  is  to  extend  anti  improve  the  "exact"  model  developeu  in 
Levine,  et.  aU,  '  in  reterence  T  to  tlie  (xiint  where  it  can  be  used  bv  motor  designers  as  a 
tool  to  aiil  in  tlie  etticient  resolution  ot  such  nonline.ir  longitudinal  instabilitv  (.ropleas. 
I  he  tirst  phase  ot  this  researcii,  tiie  results  ot  winch  are  reixirtmi  iierein,  w.is  deviileo  to 
imtirovihg  the  tiniti’  ditterenci'  numerical  technique  used  to  solve  the  equations  governiiig 
wave  [troftagat  ion  in  the  comiiustion  cliaintn'r. 

In  order  to  tx;  accefitalile  tor  the  intended  appi  ic.it  ion,  a  unite  oiilereiue  tei  hiiiniie 
must  (a)  preserve  tlit>  liigli  - 1  requeticy  content  ot  tne  w.ivelorms;  'ii,  be  ri>)ativel\ 
hondiss ipat  ive  and  non -liisfK'rsive  after  many  wave  cycles,  (i  .  be  c.rpalile  ot  desiriiang  a 
shr>ck  wave  as  a  siiarp  discont  inuitv  .vithout  ,;ener.iting  overshiKits  or  undershoots  ii(i(in 
crossing  the  discontinuity;  and  (d)  oe  cajiable  ot  p.irti.d  retlection  anu  t  r.insmission  at 

area  discontinuities.  It  should  lie  (Kiintid  out  that  in  solving  a  lombustion  instabilitv 

liroblem,  nuine  rictil  ly  induced  (in."  and  ixist-sboik  "wiggles"  oo  not  |ust  impair  the 

accuracy  ot  the  solution  but  could  also  talsidv  "trigger"  nonlinear  instabilities  aim  tone 

the  scheme  to  (lick  a  noniifivsicai  S'rlution. 

It  is  very  [litticult  tor  any  single  linite  dilterence  sclieme  to  salislv  ,ill  ot  die 

a  to  reoient  ioned  requirements  simultaneously.  Ir.r  example,  several  artitiiial  visrosiu 
schemes  havi’  been  developed  to  damp  pre-  .md  (ost-sboik  osi,  1 1  hit  imi  s.  itcnvivi'r,  sui  ii 
.1  rt  1 1  ic  ia  I  ly  introduced  diltusion  also  sme.trs  out  the  discont  imiit  s  and  eventu.ilK  iiaia|  s  the 

h  igh  - 1  requency  iilixles  that  are  [lart  ot  tlie  pliysic.ll  Ilow  lield  inside  the  motor.  Moreover, 

such  .rrtiticial  damping  i,,m  he  compar.ihle  in  .iia;,’,nit  ude  to  the  usu.il  net  gains  or  losses  ot 

acoustic  energy  in  rocket  motors,  neni  e,  its  preseni  e  would  seriousk  limit  the  validitv  ol 
the  results.  In  .uldition,  use  ot  an  .irtitici.il  vismsity  wonlo  h.imper  ettorts  to  delernane 

>.  I  ('Vine  .mil  i  iil  ii  k ,  i  ij;.  lit. 

I).  (  aplOus,  \.,  "  \  lietiiched  Stior  k  (  .1  Ir  ul.it  ion  h\  Ser  ond  I'Tder  I  mite  I  1 1 1  t  e  r  eni  es  ‘  j.  ..t 

I  omp.  Ilrys.,  kol.  .!,  pp.  Ihli,’. 


i 
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tiu‘  dctu<il  ikiinf/in^  ul  l^.is  pli.is"  in  metal  1  i/cd  Miliii  |)ri>(x‘lUmt  r<K  Let  tuotufr,. 

In  this  COP"’  ;  it  stluuld  ih-  (Joiiiti'iS  dut  that  shot  h  - 1  itt  inp  schcnit's  that  t  re.it  tHr 

shock  .IS  an  intern.il  Ijoundary  are  iinpract  ital  tor  tins  apfil  uat  ion,  (iiie  to  the  l.irye  nn.nl.er 
ot  shocks  traveling,  inte  r.ict  in^;,  and  retleclin^  insitle  variaiile  i  ros>  -  set  t  ion.i  1  <irea  motors, 
''imilarlv,  tinite  lii  t  te  rent:t‘  schemes  (such  as  the  ‘‘-scheme  or  tlie  spl  it  - 1  oe  1 1  it  imt 

i) 

.netiioti  )  that  .ire  ntKieled  alter  and  exploit  the  matliematical  tneors  ot  the  metiioo  ot 
t  ha  racte  rist  ic  s  art'  im|)ractical  tor  our  st.-ecitic  (irohle'ii.  slioik -capturing  im[ilKit 
ditterence  schemes  -itter  no  ji.i  rt  it  nia  r  aijvaiitape  Itir  tlie  iiresenl  proiile  m,  siine  the 
pinsic.il  [ir.ililem  ot  interi'st  repuires  time  resolutit)n  isinsistent  aitli  the  s[ainlit\ 
restrictions  ot  explicit  metiioi.s. 


MX  of  sesen  ee.irs  apo,  the  t.isk  >.d  I  indin;;  .1  s.iita'ole  tinite  ijillerenie  stliei.ie  .sould 

h.ive  Ik-i'ii  xirtu.ilK  iinposs liile,  .is  tlien  exist  in^  meti.ods  sut  li  .is 

Id  ,  11  Id  ,  ,  , 

'  \  .u  t  o  rm.H  k ,  l.ix  -  denur.jt  t ,  '..odunov,  *“tc.,  I.ickeo  the  rt^iuireu  pii.nities.  suite 

then,  hosses'er,  sever, il  stu'Cial  stlienu-s  tiesi^tned  t.i  .u  laeve  siiock  resolution  e.ithoul 

■  vs  i,,;files'  ssere  tlevr'loped.  Xmoiiis  tnem  .ire  the  Mux tirrei  ted- I  r.msixirt  scheme  ot  «iris 
1114 

anti  iitKiK,  Oiorin's  iintilement.it  ion  .it  fjljmm's  methotl,  tne  liistre.im  enterml  st  iieme 


7.  Moretti,  t..,  "  1  he  (  lioice  ot  .1  line  1  i.'iHMulenl  lethnii|ue  in  ;..is  Ikn.imics,"  Mh\l 

Kepirt  nn-dti,  I'olytechnic.il  Institute  ot  oriKiklyn,  Inly  Idt.d. 

t).  storetti,  "  I  h(>  );-Ssheme,'  <■  omtiuter  amt  I  liinls,  \ol.  17,  1 ')7d,  p)).  I'M -jus. 

h.  t  h.ik  r.iv.i  rt'm  ,  S.  K.,  Siuii' rsoe,  1  >.  .ind  S.ilas,  M.  I  "  llie  S|  i|  it  - 1  oe  1 1  it  ieiit  'l.itrix 

Methoii  tor  I  i\  ihtIioI  IC  Ssstems  ol  t’..is  IKnamic  I  1  |u.il  ions, '  -M  \  \  (latx'r  lUi  -  ( 1  J(  ill,  iiresent.'ii 
.It  the  llUli  Aerosii.ice  Stientes  Meelinj;,  r.is.iden.i,  t  .1 1  ilorn  i.i,  l.inu.irs  14-  In.  I'lhd. 

Id.  '.1.1 1  t  ( )  rni.ick ,  t\.  d.,  ”  I’rticei’tliiij’.s  ol  the  Settind  lnt<>ni<il  lon.il  tnntereiue  on  Numern.d 
inetnoits  in  I  IukI  1 'y  ii.im  it  s  ,  “  lertnre  notes  in  lh\sics_  ( '  t.  Molt,  1  to  ,  Sol.  ti, 

Sfir  lilt;!’ r -S  e  rlap ,  New  Sork,  I‘l7l. 

II.  Mix,  I  1.  .ind  oentlr.itl,  1..,  "Ssstem  .it  t  onse  rv.il  mn  1  .iv.  s,"  t  omm.  I'ure  Sppl.  'v'.dli., 
Sol.  I  1,  I'tnil,  p|i.  J  Ih-  J  M'. 

Id.  I  lOtitinov,  X.  h.,  "Mnite  ii|Meri'iue  Stethot.s  lor  ‘sumerit.il  '  01 1  iput  ,it  ion  s  ol 

1  '  I  St  lint  iniioiis  Solutions  ot  li|u.iti(ins  ot  I  hiiii  liMi.iimts,"  st.it'i.  s\  s.  Sol.  4  7,  I'l.'i, 

I'p.  d/l-d'ii. 

15.  !  ill  r  I  s,  I .  r.  .11  ul  I  lOok ,  II.  I  . ,  "  I  hix  (  o  1  ret  It'll  I  r.instm  1 1  111,  t\  in  imu'  n  i  r  n  1 1  I  i 

S 1  t'.or  it  hm  s  , '  1,  td  t  0111(1.  I'hysu  s,  Sol.  dd,  I't/ti,  |Hi.  5'l7-4  5I, 


14.  I  ht ir  in,  \.  | . ,  '  K,mt Ion  1 
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1 1\  rl >1  )l  u  St  "  1 .  ( » t  ;  ( )M  •.  f  1 1\  -s  u  's , 


ot  Van-Leer,*^  tlie  Hybrid  Scheitie  of  tiarteii  aiul  /vvas,^^’  and  tlie  Artificial  Coniprossion 
1  7 

■Scbenie  ot  Harten.  All  or  tliese  scfienies  arc*  finite  difference  approxiniat  ions  to  tiu* 
tierivatives  arising  in  the  conservation  lavvs  and  can  treat  an  arletrary  system  oi 
conservation  laws.  In  addition,  tliey  are  all  fixed  grid  nietfiods  and  can  automat icallv 
handle  interactions  between  waves  of  different  families.  In  orrJer  to  evaluate  the 
Suitability  of  tne  aforementioned  technirjues  for  the  present  purixrse,  their  aljility  to  treat 
a  number  of  simpler,  but  related,  problems  was  examined.  I  be  results  of  this  investigation 
are  summari/ed  in  Section  3.  I  he  finite  rfifference  scheme  |udgi*tf  most  (iromising  was 
then  incorfxjrated  into  thr?  nonlinear  instability  (rrogram  described  by  Levine  and  t.ulick^^. 
The  results  of  a  numl)er  of  nonlinear  instability  solutions  are  presentetl  to  demonstrate  ttie 
effectiveness  ot  the  new  technique. 


15.  Van  leer,  b.,  "towards  the  Litimate  (.Conservative  Hifterence  Scheme  III, 
bpstream-bentered  L  inite -Hit  ference  Scheme*  for  Ideal  Compressible  flows,"  |.  ot  Comp. 
Ptiysics,  Vol.  3,  1977,  pp.  2()l-275. 

Harten,  A.  and  /was,  (j.,  "Self  Adjusting  Hybrid  Schemes  for  Shock  (.omjiutations,"  |. 
of  Comj).  Hiysics,  Vol.  9,  1V7  2,  jiji.  5(>f5-5}ff. 

17.  Harten,  A.,  "  I  he  Artificial  t:omprt*ssion  Method  for  Conij)utat  Ion  of  Slu)cks  and  Contact 
Hiscontinuities:  III,  Self  Adjusting  Hybrid  Scliemes,"  AfOSK  lechnical  Kejxirt  lK-77-d()5d, 

March  1977. 


1(i.  Levine  and  (Culick,  Hj).  Cit. 
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I  LSI  I'KDIJLLM  KLSUl  I  S 


2.1  Sil(.x:K  IUI5L  PK11HLLM 

•\s  a  first  step  in  evaluating  the  various  tinite  tiifterente  schemes,  they  .vere  used  t<j 
solve  the  shock  tube  problem  for  the  otie -dimensional  Lulerian  form  of  the  i>as  dynamic 
conservation  equations  for  an  inviscid,  non -heat -contiuct  iny,  fluid.  Ihese  solutions  ivere 
used  to  rate  the  difference  schemes  based  u(M>n  such  criteria  as  resolution  of  the  shock 
(i.e.,  number  of  mesh  pr>ints  needed  to  describe  the  shock  iliscontinuity),  diffusion  and 
smearing  ol  the  shock  and  the  contact  discontinuities  with  time,  stability,  effect  of 
Courant  number,  anil  computation  time. 

The  results  of  the  shock  tube  test  case  demonstrated  the  superiority  of  the  recently 
1  q-2  i  24  2') 

developed  techniques  over  earlier  methods.  '  based  on  the  previously  mentioned 

criteria,  the  ftest  method  was  a  combination  ot  three  technique's.  1  he  iiasic  sclieme  of 

2  ()  27 

Lax-Wendroff  was  combined  with  a  llydrid  scfieme  and  the  Artificial  Compression 

2  (j 

Method.  fhis  methtxf  (LW +1 1  ♦ACM  J  was  capatile  of  producing  very  sharp  shtKk  isaves, 
witnout  pre-  or  txjst -shock  wiggles. 

2.2  LINLAK  WAVL  h(,)UAIION 

I  he  various  finite  difference  irteth(xis  were  also  employed  to  solve  .i  second  problem, 
i.e.,  the  linear  wave  equation.  The  results  were  used  to  evaluate  the  relative  dittusive  anu 
dispersive  errors  of  the  schemes  tor  harimNiic  st.inding  oave  protiiigat  ion  after  many  w,iv(> 


Id.  boris,  ()(>•  *  it. 

20.  Chorin,  Of).  Cit. 

21.  Van  l  eer,  Oj).  (Ot. 

22.  ifarten  and  2was,  <4).  C.it. 
2  1.  ffarten.  Of),  ('it. 

24.  Mact  ormack,  Of),  (  it. 

25.  (iOfiunov,  ()f).  (  it. 

2().  lax  and  VVendroff,  Of),  (at. 
27.  Ilarten  and  /wiis.  Of),  (  it. 
2h.  Oarten,  Op.  <  it. 
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cycles.  Ihis  proi)lein  vs.is  .ilso  used  to  assess  tfu-  <><is(‘  ol  ioipltsiieiitat  lou  ot  Ixjuiul.irv 
conditions  tor  the  tecliniques  considereii. 

llie  linear  wave  equation  solutions  demonstrated  that  tht'  ina);nitude  ot  the  ihitusive 
ami  dispersive  errors  varied  si>;ni1  icantly  trom  one  method  to  anottier.  I  he 
comtiination  im-thoci  ()ro(iuC(>d  excellent  solutions  to  this  proiilem,  with  little  error  eviilent 
alter  many  wave  cycles. 

2.^.  W  Wl  l'K(,)l>\GAl  ION  IN  A  ClOSl.D  ILHil 

1  he  tinal  test  leailiiij;  to  tlie  selectKin  ot  a  ludhod  to  he  utilized  in  the  nonliiif'ar 
instability  proj;ram  was  the  solution  r)t  the  one-diinensional,  ixrnlinear  hypt'rixilic  equations 
riesc  ribiny,  tinite  amplitiKle  wave  and  shock  propagation  in  a  closed  end  tube.  by 
establishing  ditterent  initial  comlitions,  it  was  possibk'  to  comiiare  results  obt.niu'd  with 
the  various  techni(|ues  tor  proldoms  ranging  trom  almost  line.ir  (small  initial  dist  urlianc  es  j 
to  highly  nonlinear  (large  tiisturbances)  conditions,  over  .1  large  luinbt'r  ot  wave  cycles. 
In  addition,  results  obtained  with  t!iese  techni(iues  were  coinji.irmi  lor  proldeas  such  as 
interactions  ot  ditterent  svaye  tamilies  and  the  idlect  ot  change  in  grid  size  and  (  ourant 
niinljer  ufxin  the  trcx(uency  content  ot  tile  sclu'ine. 

Son*-  ol  tile  results  ot  the  last  test  problem- - 1  inite  aniplitudr'  svavt'  in  a  closed 
tube- -will  be  iiresented  to  illustrate  sonK>  ol  the  ditlerences  bt'tween  the  techniques.  Ns 
a  result  ot  the  numerical  error  associated  with  tinite  ditterence  metinxts,  eacli  techniqur* 
at:ts  as  a  numerical  filter.  I  he  ’tiltering’  ettect  ot  eai  li  of  the  techniques  is  a  dilli-rent 
lunctitm  ot  trwjuency,  mesh  size,  et(  .  lo  enhance  tin-  ability  to  liiscern  these  dil  terences, 
the  results  ot  the  test  problems  .vere  spectrally  analyzed. 

I  igures  1,  2,  and  ^  ilemonstrate  the  ettect  ot  niinerical  tetlinique  on  the  lime 

evolution  of  acoustic  pressure  amplitude  .11  an  <md  ot  a  r  losed  tulie.  I  lu'  solution',  were 

initiated  with  a  first  liarnKznic  stamiing  wave  p<>rturbat  ion  having  an  amplitiKle  ot  20  „  ot 

the  mean  tiressure.  1  he  numerical  scln‘mes  used  in  tlie  tigures  were  the  lw+ll+\l  \t 

riiethixi  (lag.  1),  tlie  Macdormack  scheme  (I  ig.  2)  aiul  tin-  Kubin  and  biirsti'in 
2'> 

schr'ine  (Mg.  5).  I  he  evolution  of  iM)st-sh(K:k  wiggles  into  erroneous  hight>r  harmonus  is 
evident  in  the  latter  two  (  ases.  Spectral  analyses  ot  llu-  results  shown  in  Mgiires  1,  J, 
and  I  are  presented  in  figures  4,  S,  and  (>  respectively. 


2d.  Kuhin,  1.  I  .  and  llurstein,  s.  'I >it  tereiu  i-  vlethoils  tor  liu'  Invisiid  and  \  isrous 

l(|u<itions  of  a  ( onuiressihle  (i.is,"  iourn.d  ot  <  oinmtat  lonal  I’hssus,  \ol.  2,  1'hi7,  pp. 


ijjijre  2  a).  1 1'Tte  Iv.jljtiun  of  Pressure  Oscilllations  at  an  Ttid  of  a  Tube  ( MacC^onnack ) 


Unure  2(b)  and  (c). 


txpatulHd  Views  <>l  the  lime  Ivolutinn  ot  the  I’ressiin 
0s(.  illations  at  an  1  nd  id  a  I  uIm'  (Mai  (  ormat  k  ). 


jrc  i  }  .  I  I'm-  i  vitilijT  I'in  i)t  ''’ressur*'  Oscillations  at  an  l.nd  of  a  Tub**  (Rubin  and  Uurstein) 


figure  J(b)  and  (c). 


fcxpanded  Views  of  the  I  inne  I  vKiutiijn  of  Pressure 
Oscillatiijns  at  an  I.nd  of  a  Tulw'  (Kubin  and  Hurstein) 


PERCENTAGE  OF  PSD 


Figure  4(e; 
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Vurn;tii)n  of  t- renut’^^Y  and  Hurstein) 


PHRCCMTACF  OF  PSD 


^  i(4urp  h(e). 


I  irne  I  volution  of  \t  i  iimul.it  ive  I’M)  as  a  I  unit  ion  of  '.Iniln  NuuiUt 
(Kubin  and  Hurstein). 


vVilli  tile  lA+i(+\t\l  technique,  the  initi.il  line  >\,ive  <i  i  U  U  rh.ince  lieveliips  into  <i 
sharj),  t  rianpul.ir,  shock -tv;)i‘  wave  lorin  ( I  ij;.  1;.  Iliis  Is  the  expecti'ci  result  iiased  on 
the  closed  torin  .inaUtical  solution  to  this  problen,  according  to  Morse  and  lii^;ard.'^ 
Ilfiures  4d  to  4d  present  tht‘  time  variation  ot  the  /lower  s/iectral  density  contained  in 
each  .node.  I  hese  ti^ures  indicate  that  higher  inorles  u/i  to  the  l  lth  are  excited.  Ii^ure 
4e  shows  the  time  variation  ot  tlie  ac<  umulat  i  ve  (lower  s/iertral  density  anu  indicates  tliat 
tor  this  nunierical  scheiiu'  there  are  no  erroneous  shilts  ot  (mwe  r  speitral  densit\'  in  the 
dul'  order  modes.  I  he  acoustic  enerf;v  distribution  varies  little  .v  ith  time  ..iiice  .i  shock  is 
tcirined;,  .as  extiected.  1  he  transter  ot  eiu-r^v  t  roni  the  tundaniental  into  the  higher  mixles 
results  t  roll!  the  /ihysical  (irocess  ot  w.ive  stiM'/ieniii};.  Iheri'  Is  no  phy  sical  process  tor 
t  rails  te  rrin^;  enerjjv  the  other  way,  i.e,  ffoni  higher  modes  into  the  i  unmiinent  a  I. 

tijture  1  shows  the  time  evolution  ot  acoustic  pressure  ain(ilitude  at  .in  end  ot  ,i 
closeti  fuiie  '.vh(>n  utili/ing  Mact.ormack's  scheme.  \n  exjianded  view  ot  the  .icoust  ic 

(iressure  am/ilitude  lietween  nondimensional  times  ot  zero  and  titteen  is  shown  in  li^ure  dli. 
I  he  develot-iir.ent  ot  the  initial  ;xist  shock  "wiegle,"  that  at)|)ears  after  the  tirst  cycle,  into 
erroneous  hi(;her  larxies  is  evident.  -Xn  ex(ianded  view  ot  the  .icoustic  (irc’ssure  ain/ilitude 
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harm.xiics,  lias  now 

bec‘n 

reduced  to  a 

vvavetorm  c:om(iosed  ot 

4-5 

harmonics  ' 

only. 

1  he  absence  ot  the 

hijiher 

harmonics  is 

evidt.Miced  by  the  discrete  hunifis  in  the  waveform.  I  his  trend  is  also  subst  int  lated  by  the 
S(itictral  analysis  of  tins  solution  (shown  in  f  i^s.  ha  through  he).  Mytiire  hb  indicates  an 
excessively  hieh  (iprcenta(;e  ot  acoustic  ener(>y  in  the  ei^ht'i  and  ninth  harmonics.  At  the 
enii  ot  the  test  (after  M  wave  cycles),  it  is  shown  (  I  ij;.  hd)  that  tins  excessively  fn,;h 
(lercentaiie  oi  acoustic  ener(>y  has  been  t  ransterred  to  the  litth  and  sixth  harmonics, 
'lince,  ris  (ireviously  rnentlonerJ,  there  are  no  known  (ihysical  (irocesses  that  can  cause  sui  h 
a  transition  of  energy  from  a  hi/^her  inmie  to  a  lowc-r  mixie,  this  phenomenon  is  a[)|iarently 
a  result  ot  numerical  error.  f  i(;ure  he  is  anotlu'r  way  ot  demonstrating  lln‘  same  sjiuiious 
result.  Hence,  the  time  variation  in  the*  accumulative  -power  >|ie(  t  rom  is  a  result  ot  the 
dis(iersive  error  ot  the  numerical  scheme  .vliich  causes  (irc'ssure  si/;nals  to  travel  at  the 
wron^  Sjieed. 

I  he  results  obtained  ulilizint;  the  Kuliiii  and  liurstein  scheme-  (  f  i)is  f  and  •>;  a  re 
geni-rally  similar  to  the  results  out. lined  by  utili/inn  Vt.ict  ormack's  scheme  i  I  i|;s  .!  and  h  : . 

41).  Mor>e,  I’.  M.  .iiid  In/pird,  K.  V.,  I  hc*circ‘t  ic  al  Ac  lustics,  Mct.r.i.v  Hill  Hook  (  omninc  , 
Nc-w  York,  Ihtih. 

t 


\  (. onip.iri'.oii  ot  the  ri'sults  mclu.itcs  tliat  tlii"  first  (xist-stuxk  appears  after  tiie 

tiriri!  wave  eyele  eomjr.ired  to  tlie  first  wave  eyrie  (witli  sfar  Carrniai  k ).  Also,  tlie 

j)ereeiitaj4e  ot  e'ner^;v  eoiit. lined  m  tin-  t  uiula'nt'rital  larxie  is  ^ireater,  and  tlie  i-rroneously 
e(HT(4y  in  the  hi;4her  nuKles  is  soine.vhat  less  with  the  Kuhin  anti  liurstein  st  heme.  In 
tills  eoniK'ttioM,  it  shoulti  lie  mentioned  that  vi'ry  similar  results  to  those  olitaineo  h\ 
ntili/iiip  the  Istihin  and  liurstein  silienie  twre  obtained  svith  tlie  I  a\ - Weiidrof  t  seheinio 

I  he  aeeunuilat  ive  perei'ntaj;e  of  (xiwer  spectral  i/ensily  contained  in  the  res(;ett)ve 
harnionics  tor  tht'  noiuiiniensional  tinu-  tit  7.1  to  1-t.lt‘l  is  shown  in  I  ij^ure  7  tor  live  ot 
the  luinierical  schemes  consulered.  Ihe  sii|)t>riority  ot  the  l.V+ll  c  oiiiliinat  ion  sthemi-  over 

the  stiind.ird  sclunnes  is  eviiient.  1  h*‘  addition  ot  the  artificial  compression  to  tins 

coinhination  serves  to  sharpen  shock  transitions  by  exciting  hip, her  harnionics  that  were  over 
damped  by  the  addition  ot  the  llytirid  scheme.  Ihus,  coinp.iri'd  to  the  1  Vv  aiuf  IKorul 
result,  the  I  .V +1 1 +At At  solution  has  slightly  less  enery;y  in  the  tundamental  iiiixle  and 
sliythtly  more  enerey  m  tlie  ni,;her  h.irmixiics. 

All  of  the  earlier  methods,  i.e.,  kf.ictiormack,  I  ax-VVendrott,  and  Kubiii  and  liurstein 
were  used  without  addin^t  an  .irtificlal  viscosity.  the  eftect  ot  an  .irtilicial  viscosity  on 
tlie  solution  was  lnvesti);ated  usiiiy;  Hyman's'*  priulictor -corrector  scheiiK'.  Kesnlts  utili/in;; 
this  methorf  with  tsso  dittererit  artificial  viscosity  cotH I ic ient s  are  presented  in  I  ijtiires  ii 
and  4.  Ihe  aridition  of  an  artificial  viscosity  to  a  iuim<*rical  scheme  vvas  conceivc'd  as  a 
way  to  danip  post-shock  osc  il  lat  ions.  Artificial  viscosity  nxluces  (xist  shock  oscillations  at 
the  expense  of  tiie  higher  liarmonic  components  of  the  w.ivcdorm.  I  ipure  H  shows  what 
hapfiens  vvhen  a  hieli  value  ot  artificial  viscosity  is  i.'iiifdoyed  (A  e((ual  jiiity  in  Hyman's 
methcxl).  In  tins  case  the*  hiph  .irtilicial  viscosity  (ircwents  a  shock  from  ever  lormine  and 

the  deviations  trom  a  perfect  sine  wave-  <irc'  lu'ver  lar^^e.  Ihe  results  ol  the  spectral 

analysis  ot  this  solution  (Mps.  lib  and  He)  show  the  .ibsenci'  ot  higher  harmonic  content. 
Kecluc  iiif;  the  artificial  viscosity  coed f ic lent  (  iS  =l).f,  the-  lowest  valne  at  svliich  Hyman's 
method  remains  st.ililc')  yiedds  a  mne  h  steepi-r  wavedorm  (lip.  ha),  but  oni'  whose'  hipher 
harmonic  content  Is  still  Ic-ss  than  it  should  be*  (fips.  ')b  and  he).  As  time  passes,  the* 
action  ot  tiu'  artificial  viscosity  continues  to  d<im|i  pre  feri-nt  la  II  \  the-  liiplier  harmonics, 
causmp  the  solution  to  turthc-r  depeiu'ratc'. 


fl.  Hyni.in,  l.mies  M.,  "On  Itolxist  and  \cccirate  '.tethexK  tor  thc‘  C  ah  illation  ot 

I  Ompressibic-  I  luid  Mows,"  I’.irt  1  to  be  |)ubl  islic-d. 
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As  inentioiifd  earlier,  artiticial  viscosity  scIkmih'S  ivs^uld  liampcr  i‘ft()rts  t(j  dtacrniinc 
the  actual  damping  ot  gas  phase  oscillations  in  metalized  soliti  prcjpellant  s.  lo  illustrate 
the  similarity  belwerm  the  effects  ot  particles  and  artificial  viscosity,  the  previousU 
described  closed  end  tube  pniblem  was  rnoriitieri  l-y  the  aifdition  ot  2'..,  i()article  to  g.is 

weigirt  flovs  ratio)  of  S  micron  particles.  1  he  computed  results,  obtained  isith  tlie  Kuhin 
and  burstein  scheme,  are  shown  in  I  igiire  10.  Oomjjaring  figure  10  to  I  igure  1  shows  that 
the  initial  (vist- shock  wiggles  (fig.  lOt)  is  an  ttxpanded  vie.v)  did  iKit  devcdop  into  liisc  rete 
humps  in  the  .v.ivetorin.  As  expected,  the  spectral  analysis  tor  the  particle  case  dni  not 
show  the  high  errcrnt'ous  spe'Ctral  content  in  the  higfter  harmonics.  VVith  particles  present, 
the  damping  of  the  high  - 1  repuency  content  of  the  waveform  is  due  to  a  real  ph\Mc  al 

process.  Ihis  is  in  contrast  to  the  similar,  but  nonphysical,  action  ot  an  artiticial 

viscosity. 


r  rne  [volution  of  the  Norirwlized  f)sci(latorv  Cressure  at  l.i  End  of  i  '  (osnd  I  u(x‘  filled 
kvith  Particles  t  Kiibin  and  Hurstein  Scheuie). 


o 

►H 

H 

< 

PO 

H 

U4 

g 

Vi 

CO 

UA 

sx: 

a. 

o 

w 

ts] 

M 

I 

o 


Figurt?  10(b)  and  (c). 


txpanded  Views  of  Un>  Norinali/.ed  Oscillatory  Pressure 
Amplitude  at  an  fnd  of  the  C.losod  Iut».‘  Filled  with  particles 
(Kubin  and  Uurstein  Scheme). 
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MOlOK  Si<l  I  i  I  II  ilsS 

ltu“  test  prolilcii  vikitiDiis  vs('r>'  usfd  to  r.iiik  tin-  Viirioiis  nictliods  .uioroitij;  to  thr 

prt'viOLislv  riicnt  loot'll  tritt'ri.t.  It  kv.is  t  out  lutii'il  tiuit  tiif  1  ,i\ - 'tifiulr' it  I  ♦  iitloio 

\rtitH  ial  i  omp ri'ss loll  sc  ht'iiif  s  U'lclfil  tlif  Ik'sI  results,  thus,  tins  si  iieme  v..is  iiii  ■  irimraleo 

into  the  nonlinear  instahilitv  proi;rai'.  oesi  iiheil  l)\  I  evine  ano  (  nhck. 

I  he  i|uaM -one  -  diniens  loiKi  I ,  tss^i-phase,  ei|(i.il  ions  ol  i  lot  r  m  t  >  1h  soUei,  .in-  '>do,-.h  i:i 

li'^ure  11  in  lonservative  tone.  I  he  propidhint  jiroi^Tlies  and  'lolcir  ,,eonii‘t  rv 

I  \  I  innric  al  K  perforated  ^rain  ')h.7i  ni  .!  I.s  ini  lies  I014,  '  -Ai  re  t.iken  t  roM  lec  ine  .iml 
I  111  ic  K  *  '  to  tacilitate  loinpanson  imIii  the  earlier  results.  i.elori-  preseiilin,,  tin  iisnlls  ijI 
tile  motor  ^olcltlon^,  it  shoiilo  lie  mentioned  that,  lor  the  motor,  |  roiwllanl  aao  ■  a  k' r  al  ii  1 ,, 

londitlons  uidi/ed,  a  I  ine.i  r  slahilits  an,il\sts  shows  the  lunoamenlal  iiUKie  to  Ih-  niolaiile, 
while  <|II  nl  the  hl^ner  modes  are  stahl*'.  It  shoiihl  also  he  |«anted  out  that  at  die 
present  tune  velot  it\  tonplliii;  ettects  are  not  me  limed  in  the  nonlinear  inslainhtv  pro,, ram. 

(  4;nre  I  .^a  shovss  the  pressure  time  histors  at  the  head  eno  ol  Itie  inodir  lah  iilalep 

iisiiii.,  the  1  ,V  +  M  '>1  sclieine  and  an  iniidil  first  mode  (tisturh.inr  e  .nnpl  it  in  le  ol  e..;  ,,i  im,. 

mean  pressure.  Ihe  earU  hl^torv  I  roiii  t  -  d  t'l  I  >  iionoii  lensional  limi  is  sliown  ri 

espandeii  torni  in  I  i|;nie  lili.  Ihe  I'litiiillv  sinusoidal  wave  is  seen  to  .diiiost  imi  m  d  lal  e  K 
iin(Iera,i)  transit  loll  to  <1  sharp,  osi  illation  tree,  shotk-tspe  ssa\e  ti  1  r  m.  I  ipnie  1  .’1  show, 

the  pressure  histories  Imlh  at  tiie  head  enil  and  (enter  ol  die  miiloi  lot  the  lime  niteim' 

I  -  111  to  dll.  In  this  time  interv,il,  the  ssave  le.ii  his  .1  limit  nip,  a' npl  it  ni  le.  Ivaennatioii 
III  the  plia.e  ill  t  lereiu  ,'s  lietween  the  pressure  at  odiereni  loi  ilions  ihnnates  that  die 
solution  Is  neither  a  standine  w.ive,  nor  .1  Iraveliiij;  w.we,  lull  1  (o  .iiii  lal  ion  ol  •>il'i. 

Ihe  saaie  prolilen  was  soKed  nseip,  tin'  i.'ndin  ano  I'ntslein  siiiiiie  the  ntio.! 

emp  lo\  ei  1  III  Ihe  or  I,;  Ilia  I  ( iroi;  r.i  11  1 1\  I  ev  me  and  1  nl  n  k.  '  i,;iMe  s  11a,  '  ■ ,  ai  u .  '  snow  t '  a 

head  end  pressure  histor\  .niO  illustrate  l!ie  deve|o|«nenl  ol  the  initial  ,Kist  -  shoi  ^  ai  ;.'Ii 
into  an  erroneous  setoin:  mode  os<  iMatioii, 


1  J .  I  ev ii le  and  I  n I  K  k  ,  1  ip.  i  it  . 
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NORMALIZED  PRESSURE  PERTURBATION 


t  i«urc 
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(b)  .imi  fi  )•  I  Views  of  the  lime  fvoliitioii  of  the  \orm.ili/ei) 

I  )si  il  l.itorv  I’resMin  'Miifilitiide  ,it  the  Me, id  I  nd  ol  the  M 
1st  IhmiKxtK  ,  St.indiiiH  IV.ive,  '  I’  0.4  :lv\tM»  \(  M 
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N0RMA1.I7.ED  PRESSURE  PERTURBATION 


I  tu‘  nuirki'd  impruvi’Mu'nl  ahtaiDcci  .\itli  tin-  ii<  a  Mjkjti'jn  lifnir  lor  lliu 


S'  ■ \  •  I  ' 


I  <isi‘  .va  s  ,4  r.it  1 1 V  1)14. 


Ill  tile  ctliHt  ')!  initl.il  s  li  st  url)<iina'  iiiiipl  it  U(ii',  tlu*  s.iau-  pruDii’ii  a, is  .111.. s, 

USIila  tlU'  I  .1  a  1  I  ♦  SI  hvilll'  ,lllll  lllltl.ll  lirst  MIIKic  (1 1st  U  rb.lllL'l'  .Ulrpl  ltUli''S  'll  .till!  .'.'i 


ot  the  me 

.1:1  pressure. 

1  tie  iie 

■  Id  i-uct  pressure 

historic's  tor  these  c 

li  s(  -s  .1  rt  ‘  s!  K  1 1 

1  14LI  re  s  14 

.Illcl  1  1  res; 

.ec  t  ivelv  . 

III  .ill  liiree  c 

.ises,  •1.4,  II. J,  .iiid  ll.i 

1  '<  II nt  1.)  i  r I'lf ) !  It  iiu 

shossii  in 

1  1411  r.'s  1  .p 

14,  .11  id 

1  1  rc'siic'c't  1  velv , 

the  sol ul  101 1  re. 1C  he s 

1 1  u'  Sii'  nr  !  1  n  1  It  1 '  1 

.tm()l  itude  . 

It  .lb  lilt  1.  1  J 

.,1  the 

me. in  pressure. 

libs  result  w.is  ,|uitc 

‘  SI  J  r, )  r  iMn  ^  in  v  n  • . 

lit  till'  results  iis  lesiiie  .inii  (iiliek  vvliiell  rlemiiiist  rated  <1  st  riiii^  etteit  it  initial 
cl  ist  Liriainee  iin  limit  in;.,  aniplitiKie.  Ihe  explan. it  ion  .i|ipears  to  lie  in  t  ne  tail  tli.it  ttie 
(iresent  solutions  v\ere  lor  4.1s  only,  >\hile  tlie  e.irlier  solutions  .ser.'  lor  the  s.mie  a  il'ii, 
0()er.itin4  Conditions,  .ind  p  n  ipell.int ,  but  Aitb  1 'i  ,  .ihiiununi  ad.ied  nl.bi  p.irtuh'  4.1s  .\ei„iil 
t  lov\  ratio).  llent'e,  it  seems  that  nonline. ir  p.irticle  d.iinimi^  eliecis  |ila\  a  1  nt  u  .d  role 

in  estdbl  ishiny;  limiting  anipl  it  cKle. 

Ihe  spectr.tl  .iii.tlvtiib  results  cor res|K)ndiii4  to  the  il.idi  .ii.iplitude  i  .ise  am  shown  n 
I  14Li.es  Ilia  through  Ino.  llu“se  results  <U‘.irl\  demonst  r.it  e  tiie  excitation  ot  iiiaher 
h.irnionics  .is  the  origin, tl  .Ji sturhance  jjro.vs  .ind  steepens. 

•\n  import, int  .ip(il  icat  !■  in  ot  nonllne.u  inst.ibilitv  .in.dvsis  1,  the  prediction  ot  mot.ir 

response  to  pulse  type  disturb. inces.  I  be  I  evine  .ind  ciilick^'  nonlme.ir  inst  .tin  I  it  s  tiio^ran 
h.is  .1  liiiiti'd  .iliil'tv  to  ctein'iate  pulse  tyfie  cvavelorms  it  will  lie  exii.inded  soon;.  Ibis 

cap.ibilitv  w.is  used  to  test  the  .ibility  ot  the  |  .vs||+\(  \\  method  to  tre.it  such  problen-. 
1  lie  results  ot  t  .VO  solutions  .vill  pe  discussed.  In  imth  cases,  the  initial  iiri'ssure 

disturb. nice  .v.ivelorm  was  taken  to  iie  ot  the  lorm  sin  1  ^  vvith  ,in  amplitude  ei)U.d  to 

U.4  ot  the  me. in  pressur.'.  I  his  pr'iduces  the  centered,  symincM  ric  evavetorm  sliocvn  in 

I  inure  1'.  Ihe  ditlereium  iiet-veen  tite  t -\o  c.ises  was  the  initial  velocitv  .11  t  -  l).  In 

one  (  .ise,  tile  veloc  its  A. is  t.iken  to  be  !’/  .  vvhile  in  the  second  c.ise,  the  velocitv 

V  nondii. lens  i.  III. d  ...is  t.iken  t.i  be  /ero.  Ilte  lirst  c.ise  rc'presents  .i  travelin;;  pulse.  I  lie 

vahie  ot  veloc  itv  is  i.ikeii  sue  h  th.it  the  initial  (lulse  represents  ,1  rislit  tr.lvelinn  A.ive.* 
Ihe  second  c.ise  c  or  ri'S[)oiids  to  .1  "  st.i  lid  illy;’  pulse'.  Ihe  pulse  p  rot). 14.it  es  .is  the'  sun  cd 

epu.il  Ic'tt  .iiid  r  lyht  t  r.ivc-l  1114  vv.ivc's,  e.icii  haviu4  halt  the  initi.d  .im|il  itude. 

d>.  Iind. 
bi.  loid. 

•■'sctii.illy  settle^;  V  0  '  p,  ,  only  produces  ,1  pure'  ,iud  riplit  tr.ivcdiui;  .s.ive  in  I  ne  line.ir 

limit  IS  '  I’  .ipiiro.ic  lies  1). 


9*6 


I hf  cdlculatecl  |)ri‘ssi.ir<‘  lif>t()r  »‘s  dt  the  lu'di)  iMui  nt  the  I  n  mK  ii  .n  iIum' 

ilisturbdiu  ns  dm  shown  in  t  r^iirc  '  ;i  Ufdvelinj;l  dii<l  I  i^uri'  in  sl.inoinn  .  Mn-  iirin.ith 

di  t  leronti'  hi'twin'ii  the  ri-siilts  ilfinonsl  rates  tlu-  nntiort.iiii  c  ol  tin-  sclcx  its  1. 1  n  1 1  i.mi  i  ■ 
dssocidttnl  with  a  prs'ssiire  iialse.  I  In-  travr-liiij;  pulse  is  imiiieilidl  i-K  ti.iiist  riiieo  inli 
St  ei'() "  t  font  ed  shoik  tv|)e  ssavetorii  <iiid  i|i‘<  .i\  s  until  it  ttnu  hes  the  same  liniitunt  tniil  it  ide 

as  tlie  solutions  started  t  roin  tir-.t  iiiiKie  sinusoidal  dist  uroaiu  es.  Speitral  anal,  sis  .it  inis 

solution  tor  three  t  iuu'  intervals  are  slxissi  in  I  ij;ures  iOa  to  ^Oi  .  \t  e.irls  tunes  I  i^;. 
2l)<i)  tlte  traveling  pulse  disturhaiKC  eont. lined  i  laree  |  le  ri  ent  a);e  ot  the  t  undaiiieiit  a  I,  mt  a 
siHitilitant  hl>ther  harinmic  (ontent  is  also  evident.  \s  the  limit  in,;  amplitude  was 
a()proac  t>'d,  at  later  tiiai's,  tiu'  pereent.ine  ol  in  tiu  tundamenlal  iiu  r.  aseo. 

I  In-  prt'Ssure  history  ot  the  standiiij;  pulse’  de.turhaiu  e  is  shown  in  lis;ure  In.  itie 
tune  variaiion  i  I  the  svavetorin  is  unite’  eoin(ile’V  in  this  e  .ise-.  I  lie  sjn’eti.il  mah,  sis 

results  sheiwn  in  I  i^ure’s  i  la  to  k  he’lp  to  clarity  wli.it  is  happe’inny.  \t  e.irk  time  s,  a 
syrnitK’tric  standiii).',  pulse  ceiite*reel  in  the’  motor  eontams  esseiitialK  onl\  even  h.irnione  s 
(,  I  ij',.  2  1a),  with  the  ss’cond  h.iriiuuiK  iloininat in^.  llie  1undanient.il  .ind  odd  harmimii  s  ,ire 

in f III itesiiui I  at  this  time’.  Since’  only  the  luiKiame’iital  is  uiist.ihU’  tor  this  motor,  the’  ese  i 
harmonics  dee:a\  with  time,  while  the’  linelamenl.il  hej;ins  to  i;row.  In  tiii’  time  inters  d 
f  roiti  aliout  t  =  2d  to  40,  he’  ssave’torm  hi’e.'ome’s  com|ile\  .is  it  transitions  tn.'ii  ,i  steep 
si’e  eiiiel  doiiiinati’il  vv.ivi’  to  an  .llniost  sinuso’eJal  w.iye  .it  the  t  uiid.in’ent  .1 1  t  ri’i|ue’ni  1  . 
Ke’cently,  this  solution  w.is  eoelinue’el  eieit  to  .1  noiielinii’iisioii.il  time  ot  liVi.  1  he’  ss.ise 
eonlinui’s  to  eleeiay  out  to  .1  time’  ol  alioul  OH).  \t  this  time’,  the’  w.ive  .implitude  is  oni, 
l.ti,,  ot  the’  me’. in  (iressuri’  (com|>.ire’ei  to  40.,  mitl.illv)  .inel  the  w.ivi’  is  I’ssi’iit  la  II  \  .1  pure 

t unel.i'iient .il  sine'  svave’.  \tter  t  —  100,  the  w.ive’  st.irts  to  ;;rovs,  .iiiel  is  still  Krowini;  .0 

t  =  I  do.  \s  the’  w.ive’  itro.vs,  hi,4lM’r  h.irmonic  eonte’iit  .i).;.iin  lie^ins  to  .ippe.ir  ,is  .1  result 

ol  erieritv  tr.inste’r  from  the’  tunil.imenl.il.  il  is  e\p<’r.i’il  Ih.it  the-  oseill.ition  will  e  out  nine 
to  >;row  eiritil  It  ri’.ie  he’s  the’  s.ime’  limit  ill);  .inipl  it  mli’  .is  the  othi’r  test  i  .ises  tor  this  mot.ii 
.mil  pn  ipel  hint . 
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